Non-equilibrium triplet blockade in parallel coupled quantum dots 
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It is theoretically demonstrated that parallel weakly tunnel coupled quantum dots exhibit non- 
equilibrium blockade regimes caused by a full occupation in the spin triplet state, in analogy to the 
Pauli spin blockade in serially weakly coupled quantum dots. Charge tends to accumulate in the 
two-electron triplet for bias voltages that support transitions between the singlet and three-electron 
states. 
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From fundamental aspects of spin and charge corre- 
lations the two-level system in a double quantum dot 
(DQD) has recently become highly attractive. It has 
been demonstrated that spin correlations lead to Pauli 
spin blockade in serially coupled quantum dots (QDs), 
where the current is suppressed because of spin triplet 
correlations, 013, S 0] something which may be appHed 
in spin-qubit readout technologies. [i| Pauli spin block- 
ade has also been reported for general DQDs with more 
than two electrons. |6| Recently, the Pauli spin blockade 
with nearly absent singlet-triplet splitting has been em- 
ployed in studies of hyperfine couplings between electron 
and nuclear spins. d, 0,0,0,^1 Besides being present 
in serially coupled QDs, it is relevant to ask whether an 
analog of the Pauli spin blockade is obtainable in parallel 
QDs. 

The purpose of this paper is to demonstrate that par- 
allel coupled QDs, see Fig exhibit regimes of non- 
equilibrium triplet blockade. Here only one of the QDs 
is tunnel coupled to the external leads while the second 
QD functions as a perturbation to the first QD. Impor- 
tant quantities in order to find the non-equilibrium triplet 
blockade regime is that the QDs are coupled through 
charge interactions, e.g. interdot Coulomb repulsion and 
exchange interaction, and weakly through tunnelling. In 
absence of interdot exchange interaction there may be 
regimes of usual Coulomb blockade in a finite bias volt- 
age range around equilibrium. 

In presence of a sufficiently large ferromagnetic 
interdot exchange interaction the triplet states 
|(7)a|(t)b, ct =T)i (one electron in each QD with 
equal spins) and [| t)^! Vib + | 1)^1 T)b]/a/2 acquire a 
lower energy than the lowest two-electron singlet (the 
singlet states being superpositions of the Fock states 

{[I T)a| i)B- \ i)A\ T)b]/V2,| n)A\Q)BMA\n)B}). 

Then, the triplet naturally becomes the equilibrium 
ground state with a unit occupation probability, pro- 
vided that the two-electron triplet state has a lower 
energy than all other states. The triplet persists in 
being fully occupied for bias voltages smaller than the 
energy separation between the triplet and singlet states, 



although transitions between the one-electron states 
and the singlets may open for conduction. However, 
for larger bias voltages this low bias triplet blockade 
is lifted as the transitions between the triplet and the 
one-electron states become resonant with the lower of 
the chemical potentials of the leads. At this lifting, the 
current through the system is mediated via transitions 
between the two-electron singlets and the one-electron 
states. 

The non-equilibrium triplet blockade regime is entered 
at bias voltages such that transitions between three- 
electron states and, at least, one of the singlet states be- 
come resonant, see Fig. Q] while transitions between the 
three-electron states and the triplet lie out of resonance. 
At those conditions, an electron can enter the DQD from 
the lead with the higher chemical potential, through tran- 
sitions between the singlet and the three-electron states. 
Transitions from the triplet to the three-electron state are 
suppressed since the bias is lower than the energy bar- 
riers between those states. However, electron tunnelling 
from the three-electron states in the DQD to the lead 
with the lower chemical potential are supported through 
transitions from those states to the triplet, since the tun- 
nelling barrier to this lead is sufficiently low. Thereto, the 
probability for such transitions are about unity whereas 
the probability for transitions between the three-electron 
states and the singlet is at most a half. Finally, charge 



t, U', J 




FIG. 1: (Colour online) Left panel: The coupled QDs of which 
only one is tunnel coupled to the leads. Right panel: Pro- 
cesses leading to the non-equilibrium triplet blockade. Faint 
and bold lines signify low and high transition probabilities, 
respectively. See text for notation. 
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that end up in the triplet through this process is trapped 
in this state because of the neghgible probabihty for tran- 
sitions between the triplet and the one-electron states. 

It is noticed that, a finite (ferromagnetic) interdot ex- 
change interaction is not a necessary condition for the 
existence of a non-equilibrium triplet blockade regime. 
Nevertheless, a ferromagnetic exchange yields a larger 
degree of freedom in the variation of the interdot tun- 
nelling and, also, allows a higher temperature. 

For quantitative purposes, consider two single level 
QDs {ea, £b, spin-degenerate) with intradot charging en- 
ergies {Ua,Ub), which are coupled by interdot charging 
([/'), exchange (J > 0), and tunnel ling (t) interactions. 
Specifically, the DQD is modelled by[il[ll[i3 Hdqd = 

Y.^=AAT.a^^Aad^^ + Un.^Tl,^) + (U' - J/2){nA^ + 

n-Ai ) {nBi +nBi)-'2JsA-SB + J2a (td-AcrdBa + H.C.), where 
= (1/2)I]ct<t' dl^&aa'dta', cr,cr' =T,i, A,B, are the 
spins of the two levels. In analogy with the Pauli spin 
blockade in serially coupled QDs,[l|,0| it is required that 
the lowest one-electron states, the triplet, and the two 
lowest singlets are nearly aligned, and that the lowest 
three-electron states lie below the equilibrium chemical 
potential /Lt. Hence, Et ~ Esi « Es2 ~ niin^^j{i?i„} < 
iiiiii^=i{^3n} < /i < Ei, where Et and Esn, 1,2, are 
eigenenergies for the triplet and the two lowest singlet 
states, respectively, whereas Ein, E^n, and E4, are the 
energies for the one-, three-, and four-electron states, re- 
spectively. This requires that fi — eb ^ Ae, U' « Ae, 
and Ua ~ 2Ae < Ub, where Ae = eb ~ £a- The in- 
equality Ua < Ub points out that the QDs do not have 
to be identical, merely that the charging energy of the 
second QD should be lower bounded by the charging en- 
ergy of the first. It should be emphasized, however, that 



the presence of the second QD is essential in order to 
obtain the effect discussed in this paper. Finally, weakly 
coupled QDs, e.g. ^ — 2i/Ae <^ 1 implies that the ener- 
gies for the lowest one- and three-electron states acquire 
their main weight on Q,T)a- This condition yields a low 
(large) probability for transitions between the triplet and 
the lowest one-electron (three-electron) states. 

In general there are 16 eigenstates of Hdqd, labelled 
{\N,n),E]\in} denoting the nth state of the A''-electron 
{N — 1, . . . , 4) configuration at the energy i?Ar„.0| In di- 
agonal form, the DQD is thus described by Hdqd = 
^^^ENn\N,n){N,n\. Taking the leads to be free- 
electron like metals and the (single-electron) tunnelling 
between the DQD with rate Vka-, the full system can be 
written as0 

H^ ^ Ekacl^Cka +HdQD (1) 
fe(TGL,-R 

+ HAdAa)TN+AjN,n){N + l,n'\+H.c.], 

ka.Nnn' 

where (rfA(T)ArW+i ~ {N.,n\dAcr\N + f,n') is the matrix 
element for the transitions \N,n){N + l,n'|. The opera- 
tor dAa signify that electrons tunnel from molecular like 
orbitals in the DQD through QD^ to the leads, which ap- 
propriately describes the physical tunnelling processes. 

Following the procedure in Ref. \^ the occupation of 
the eigenstates are described by a density matrix p = 
{\N,n){N,n\}. In the Markovian approximation (suffi- 
cient for stationary processes) one thus derives that the 
equations for P/vn = {\N,n){N,n\) to the first order in 
the couplings T^^^' = 2Tr J^keL/R M^H'^Ska) = To/2 
between the DQD and the leads, can be written as 



^-FjVn — X/ (y^.^N-ln' ,N7i[fa {^Nn,N-ln')PN-ln' - fa {^Nn,N-ln')PNr. 



a—L.R n 



y^r^^ jV^l„, [f;^{AN + l7i',N7i)PN7i - fa (AAr+l„',Ar„)^Ar+l«'] ) = 0, 



(2) 



iV = 1, 



,4, 



where P_i„ = P^n = 0. Here, AN+in'.Nn = En+iu' - 
Ej^n denote the energies for the transitions \N,n){N + 
l,n'|, while r^^/f^+i„, = E.^'^^^'idAarNN+i- Also, 
f^^j^iuj) = /(w — Hl/r) is the Fermi function at the 
chemical potential Hl/h of the left /right (L/R) lead, and 
f^^j^^oj) = 1 — f^^jjiuj). Effects from off-diagonal oc- 
cupation numbers {\N,n){N^n'\), which only appear in 
the second order (and higher) in the couplings, are ne- 
glected since these include off-diagonal transition matrix 



elements to the second order (or higher) which generally 
are small for ^ ^ 1. 

Since the low bias triplet blockade can be found for 
weakly coupled QDs whenever J > is sufficiently large, 
the following derivation focus on the non-equilibrium 
blockade. The non-equilibrium blockade discussed here, 
is driven by opening transitions between the two- and 
three-electron states. For simplicity, assume that the 
bias voltage V = {fiL — A'-r)/^ is applied such that 
Hl/h = /I ± eV/2. Then for \eV\ < TAe/A, kBT < 
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O.OIUa, and ^ < 0.2, which is sufficient for the present 
purposes, only the population numbers Pin, n = 1,2^ 
Nt - -P2„/3, n - 1,2,3, F24, P25, and Psn, n = 1,2, 
are non-negligible. The other populations are negligi- 
ble since the corresponding transition energies lie out 
of resonance. Because of spin-degeneracy it is noted 
that Pin = Ni/2, n = 1,2, and P^n = N3/2, n = 
1,2, which reduces the system to five equations for the 
population numbers. As discussed in the introduction, 
the non-equilibrium blockade arises when transitions be- 
tween a singlet and the three-electrons state are reso- 
nant. Therefore, the bias voltage is tuned into the regime 
where lies around these transition energies, e.er.llfij 
min„„/{A3„/_2«} < t^-L < niax„„/{A3„/^2n}, n= 1,...,5, 
n' — 1,2 (here eV > 0, the case eV < follows by 
symmetry of the system). For such biases it is clear 
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/fl(A2,i,ln') ^ 1, n = 1,. 



that /^(A2„, 

n' = 1,2, and that /j^(A3„'^2ri) = 0, n = 1,...,5, 
n' = 1, 2. It is also clear that the charge accumulation 
in the triplet is lifted for biases that supports transitions 
from the triplet to the three-electron states, hence, the 
bias voltage has to be such that f^^j^{^3n',2n) ~ 0, that 
is A3„/^2n = Esn'-Er > tJ-L + ksT, n = 1,2,3, n' = 1,2. 
Thus, the equations for the population numbers can be 
written as 



Ni 

P2n 
P 



2/3 



Nt 



(3a) 



p ^ \ + 2p{nlPY 
l^i±A|pE^^(A3i^ (3b) 
2 L2+A2/+(A3i,2„) ' ' ^ ^ 

A2/+(A31,2„ 



n— 4 



Ll + Klft{^,lan) 



3k^ 



Here, (n' 

em + 



x/„ (A3i,2„) 



1,2, n 



1 - 



A,i/^(A31,2„) 
L2+A2/+(A31,2„) 



(3c) 



\n'n\2 ^2 



4,5) 13^ 

l(rfA.)^3"'P = (1 + e)/[{l + + ei and 

A^ = J2a\i'^Aa)23 P are the matrix elements for the 
relevant transitions. The above relations are due to spin- 
degeneracy, e.g. A2„,ll = A2„,12 and A3i,2n = A32,2n, 

n = 1,...,5. Using Eq. charge conservation (1 = 
T,Nn Pnu =Ni+Nt + E„ P2n + N3) thus imphcs that 



N^ 



1 



2/3 



1 + 2p{K/py 



l+p 



l^ Ll+pAlZ^f-jA^i^.n) 

L2+A2/+(A31,2„) 



2 . 

n=4 



(4) 



Now, the matrix elements L^,A^, n = 4,5, are finite 
and bounded, however, L4,2A5 1 and Ll,Al 
as ^ — ^ 0, hence, the last term in Eq. 10} is at most 
1/2 for weakly coupled QDs since p ^ 0, ^ ^ 0, in the 
considered bias regime (see discussion below). However, 
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FIG. 2: (Colour online) Variation of tlie ratio 2p{k//J)^ as 



function of J for different t at constant Ae, U' , Ua/b- 

^2 



The 



inset shows the region in {t, J)-space where 2p{hi/ (5) > 10 



the ratio 2p(k//3)^ is finite for all ^ and J > 0, while it 
diverges as ^ ^ for J = 0, see main panel in Fig. |21 For 
weakly coupled QDs one thus finds that Nt ~ 1/(1-1- [H- 
2p(k//3)2]~1) « 1 whenever 2p{k/(3)'^ > 1. The inset of 
Fig. 1^1 illustrates a subset in {t, J)-space where this ratio 
is larger than 10^. At this condition, the boundary is 
approximately given by J{t) = Jq — 15t^[l -I- (lOi)^]. 

Using the transport equation derived in Ref. 17, iden- 
tifying G^„jY^;^„,(w) = i2TrPN+in'S{uj-AN+in'.Nn) and 

GNn,N+in'i^) = -i^TtPNnSiuj - An +in' ,Nn) , the Current 
in the considered regime is given by 



I = 



6h 



3{(3'-^') + J2[Ll-AlfZiA3i,2n)] 



_L0V^A2f+^A A + P^l J2a fa (A31,2n) 
+ ^ 2^ Anh (A31,2n) + ^ 



Nt 



1 4- 2p{>i/py 



(5) 



This expression clearly shows that a large value of 
2p{k/P)'^ yields a suppression of the current, that is, at 
the formation of a unit occupation in the triplet state. 
For biases such that fi^ < min„„/{A3„^2n} is follows that 
/^(A3„/^2n) ~ p ~ 0, which accounts for a lifting 
of the triplet blockade where the current is ~ 2p(k//3)^ 
larger than in the blockaded regime. 

The non-equilibrium triplet blockade depends on the 
interplay between J and t. A reduced t leads to a strong 
localisation of the odd number states in either of the 
QDs, which for Ae > leads to that the lowest odd 
number states are strongly localised on QD^. Then, the 
probability for transitions between the triplet, and the 
one-/three-electron states is small/large (/? ^ 0/k ^ 1). 

The singlets, on the other hand, are expanded in terms 
of the Fock states {[| T)a| 1)b - \ 1}a\ T)b]/V2, I Ti 
)a|0)b, |0)yi| Ti)-B} with weights that are slowly vary- 
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FIG. 3: (Colour online) Variation of the triplet occupation 
number A^t a) and the modulus of the current (units of eFo / h) 
b) as function of the bias voltage and the equilibrium chemical 
potential jj,. Here, ^ = 0.01, fcsT = O.OIUa = 4t, and J = 
0.2{Ua - U')/2. 



ing functions of t, however, strongly dependent on J. A 
negligible J yields that the two lowest singlets are al- 
most equally weighted on the states [| t)^! Db — \ Da] T 
)b]/V2 and | ti)A|0)_B. Increasing J > redistributes 
the weights such that the lowest singlet (|2,4)) acquires 
an increasing weight on | ti)A|0)_B, whereas the second 
singlet (|2,5)) becomes stronger weighted on [| 1)^1 i 
)s - U)a| T) b]/V2. Hence, for a finite J > and t 0, 
this redistribution leads to that transitions between the 
lowest one-electron states and |2, 4) (|2, 5)) occur with an 
enhanced (reduced) probability, e.g. L| —> 1, (L| 0), 
and oppositely for transitions between the singlets and 
the three-electron states, e.g. A4 0, (Ag 1/2). This 
implies that p as t — > while p{k/(3)'^ remains al- 
most constant. This constant, however, becomes larger 
(smaller) for smaller (larger) J. 

The typical variation of the triplet state occupation 
number iVr, calculated from Eq. as function of the 
bias voltage and the equilibrium chemical potential for 
< J < Jo - 15i2[l + (10t)2] and t/ikeT) < 2 is plot- 
ted in Fig. 131 a). Here, varying the equilibrium chemical 
potential mimics the effect of applying an external gate 
voltage Vg by means of which the levels of the DQD are 
shifted relatively the equilibrium chemical potential. The 
extended diamond marks the region where the occupa- 
tion of the triplet is nearly unity and where the transport 
through the DQD is blockaded. The calculated current 
is displayed in Fig. O b) , from which it is legible that 
the triplet blockade regime is subset of a larger domain 
of a nearly vanishing current through the DQD. The two 
diamonds within the low current regime are caused by 
a lifting of the triplet blockade (see the introduction), 
where the current is mediated by transitions between the 
one-electron states and the singlets. 



As is seen in Fig. |3| shifting fi in the range Eb + 
(Ae — J, 2Ae) causes an extension of the low bias triplet 
regime since the transitions between the triplet and the 
one-electron states become resonant at higher biases. On 
the other hand, the non-equilibrium triplet blockade is 
shifted to lower biases since n lies closer to the transition 
between the three-electron states and the singlets. The 
two blockade regimes merge into single as \fj, — A3n',2n\ < 
1^ - A2i,i„'|i n = 4,5, n' = 1,2, e.g. for fi - eb ^ 
(3Ae/2, 2A£), see Fig. O Shifting /i in the interval es -I- 
(Ae/2, Ae — J) removes the low bias blockade since the 
one-particle states become the equilibrium ground state. 
The non-equilibrium blockade is shifted to lower biases, 
here caused by transitions between the one- and two- 
electron states which tend to accumulate the occupation 
in the triplet. 

While the case Ae > is considered here, the non- 
equilibrium blockade is also found in the opposite case, 
e.g. Ae < and fi — ea ~ Ae. In this case, however, the 
system has to be gated such that only the four-electron 
state lies above /z, whereas the charge accumulation of 
the triplet state is governed by the same processes as 
described here. 

It should be noted that higher order effects, as well 
as singlet-triplet relaxation, have been neglected in the 
equation for the population probabilities Pnu- However, 
in many aspects the situation discussed here corresponds 
to the experiment reported in Ref. Q, hence the effect 
considered should be measurable under much the same 
conditions. Therefore, as in the case of the serially cou- 
pled DQD, the higher order effects give contributions 
that are at least two orders of magnitude smaller than 
the second order contributions. Therefore, these can be 
neglected in the present study. On the same basis as 
in the description of the serially coupled DQD,0 the 
singlet-triplet relaxation may be neglected here. 

The conditions required for the existence of non- 
equilibrium triplet blockade, concerning the intra- and 
interdot charge interactions for weakly coupled QDs, 
have been experimentally obtained for serially coupled 
QDs.0, 0, m The additional requirement, i.e. a ferro- 
magnetic interdot exchange interaction which is larger 
than the interdot tunnelling and the thermal excitation 
energy, is accessible within the present state-of-the-art 
technology, j3ii. il 2nl8i. il a] 
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